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A Brownian motion model is devised on the manifold of 
S-matrices, and applied to the calculation of conductance- 
conductance correlations and of the weak localization peak. 
The model predicts that (i) the correlation function in B has 
the same shape and width as the weak localization peak; (ii) 
the functions behave as oc 1 — 0(B 2 ), thus excluding a linear 
line shape; and (iii) their width increases as the square root of 
the number of channels in the leads. Some of these predictions 
agree with experiment and with other calculations only in the 
limit of small B and a large number of channels. 
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I. INTRODUCTION 

The experimental observation of quantum interference 
effects in submicron devices has triggered rapidly growa 
ing interest in the properties of mesoscopic systems.u 
Many investigations have focused on electronic trans- 
port through such microstructures and, in particular, on 
the behavior of the conductance G(X) as a function of 
some parameter X. (Two examples for X are the Fermi 
wave vector k or an external, time-independent magnetic 
field B.) Both in diffusive transport through quasi-one- 
dimensional disordered structuresa and (more_recently) 
in ballistic transport through chaotic cavitiesn one has 
observed the following two phenomena, (i) As the pa- 
rameter is being varied, the conductance fluctuates ir- 
regularly around its mean value. These fluctuations have 
a typical magnitude of the order e 2 /h. This magnitude 
is independent of material, geometry or degree of disor- 
der of the probe ('universal conductance fluctuations'); 
however, it does change depending upon the presence 
or absence of time-reversal symmetry, (ii) The average 
conductance at zero magnetic field is consistently lower 
than the average conductance at large magnetic field, 
the difference being again of the order e 2 /h ('weak lo- 
calization'). These two phenomena- are well accounted 
for not only by microscopic the£>cies,l3 but also by macro- 
scopic random matrix theories Jj~Q The success of the lat- 
ter clearly demonstrates the universal nature of the ob- 
served effects. 

The universality of the magnitude of conductance fluc- 
tuations and of the weak localization correction being es- 
tablished, the obvious question arises whether also the 



frequency spectrum ('power spectrum') of the conduc- 
tance fluctuations, as well as the width and shape of 
the weak localization peak, are universal. The power 
spectrum is the Fourier transform of the conductance- 
conductance correlation function 



C(AX) := (5G(X)6G(X + AX)) 



(1) 



where 5G := G — (G) denotes the deviation of the con- 
ductance from its mean; while width and shape of the 
weak localization peak are encoded in the function 



W(B) := (G)(B)-(G)(B = oc) 



(2) 



(The averages are taken over an ensemble of probes with 
in the case of disordered structures - different impu- 
rity configurations, or - in the case of chaotic cavities - 
different shapes. By a widely accepted ergodicity hypoth- 
esis, ensemble averages are equivalent to the average over 
some parameter, wch as k or B.) Both C(0) and W(0) 
are well knownpta the objects to be investigated are the 
ratios C(AX)/C(0) and W(B)/W(0). 

In this paper I focus on the ballistic regime; moreover, I 
assume the ballistic cavities to be classically chaotic. For 
such systems the ratios C{AX)/C{Q) and W(B)/W(0) 
have been calculated semiclassically. For example, semi- 
classical analysis yields a Lorentzian, 



C(Ak)/C(0) - [1 + (Afc/ 7cl ) 2 



(3) 



or a Lorentzian-squared, 



C(AB)/C(0) = [l + (AB/aM 2 }- 2 , (4) 

respectively, for the correlation function!! (Here 7 c i and 
a c \ are characteristic inverse length and characteristic in- 
verse enclosed area, respectively, of the classical trajec- 
tories; 0o = hc/e is the elementary flux quantum.) The 
semiclassical result for C(AB)/C(Q) has been checked 
experimentally, and appears to be in excellent agreement 
with the dataH For the weak localization peak, semiclas- 
sical theory predicts a Lorentzian,Ej 



W(B)/W(0) = [l + 2(B/a 



(5) 



This prediction, too, agree^s/ell with experiments and 
with numerical simulations nu We observe that the cor- 
relation function in B and the weak localization peak 
look almost identical. Indeed, for small B 



C(B)/C(0) « W(B)/W(0) « 1 - 2(B/a cl fo) 



(6) 
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only the tails at large B differ. 

The above semiclassical calculations involve various 
approximations. In addition to (i) the semiclassical limit 
proper, it is assumed that (ii) the magnetic field is suf- 
ficiently weak so as to affect only phases, but not the 
classical trajectories, and that (iii) in the case of weak lo- 
calization, only symmetry-related paths contribute ('di- 
agonal approximation'). The first two assumptions are 
expected to be justified only in the limit of weak fields 
and a large number of channels; while the third assump- 
tion has in faet been shown to be too crude.El Moreover, a 
recent studyli3 has revealed that important contributions 
to weak localization stem from diagrams which appear to 
have no semiclassical analogues. In view of these theoret- 
ical limitations it is quite remarkable that semiclassics is 
so successful in describing the shapes of C and W . Nev- 
ertheless, it seems highly desirable to derive parametric 
correlation functions and the weak localization peak by 
some alternative method. 

The above results (||), (Q) and (||) are universal in 
the sense that the geometry of the probe enters only 
via the characteristic scale 7 c i or a c \. Such universal- 
ity suggests that the phenomena in question should lend 
themselves to a treatment in the fraMework of random 
matrix theory. Indeed, Pluhar et al.til showed that the 
Lorentzian shape of the weak localization peak can be 
derived with Hamiltonian random matrix theory. In ad- 
dition, they succeeded in proving that the width of the 
weak localization peak increases essentially as the square 
root of the number of channels. (This is a result which 
in a semiclassical theory can be inferred onlv_indirectly.) 
As regards parametric correlations, AltlandE3 employed 
Hamiltonian random matrix theory to show that in the 
diffusive regime and quasi-one-dimensional limit, C(AE) 
is a Lorentzian. (E denotes the Fermi energy of the elec- 
trons.) Other macroscopic calculations of parametric cor- 
relations have so far focused on correlations of the level 
density, rather than on correlations of the conductance.tS 
Their success certainly encourages further applications of 
macroscopic random matrix models. 

A particularly simple and (arguably) generic model 
for any kind of parametric dependence is furnished by 
Dyson's Brownian motion modelo It is based on the as- 
sumption that variation of the parameter X amounts to 
a random walk on the manifold of Hamiltonians, a pro- 
cess governed by a diffusion (or Fokker-Planck) equation. 
Beenakkenla successfully employed this model to calcu- 
late density-density correlations, and suggested that it 
might also be applied to the response of transmission 
eigenvalues (and hence of the conductance) to an exter- 
nal perturbation. This is difficult, however, because the 
conductance - in contrast to the level density - cannot be 
easily expressed as a function on the manifold of Hamil- 
tonians. This difficulty can be avoided if one considers 
random S-matrices rathec than random Hamiltonians. 
Random S-matrix theory,Q while consistent with Hamil- 
tonian calculationsu in the limit of a large number of 
channels, allows a much more direct calculation of con- 



ductance properties. Therefore, I propose that the Brow- 
nian motion model be reincarnated as a random walk on 
the manifold of S'-matrices J13 To devise such a Brownian 
motion model for S'-matrices, and to study its implica- 
tions for conductance- conductance correlations and the 
weak localization peak, is the purpose of this paper. 



II. THEORETICAL PRELIMINARIES 

I will first review some basic definitions. Electronic 
transport through a microstructure constitutes a quan- 
tum mechanical scattering problem. O Scattering at a 
probe with two leads, each with TV channels, is described 
by the S-matrix 



r t' 
t r' 



(7) 



where t, t' and r, r' are the N x N transmission and re- 
flection matrices, respectively. Both (t't) and (t'H 1 ) have 
the same set of eigenvalues {r a £ [0, 1]}. In terms of these 
eigenvalues, the S-matrix may be parametrized 



S = 



with 



«w 
o «w 









(8) 



N 



=. diag(r a ) and the {»'''} unitary N x 
matrices.ta The {r a } and will be referred to as 

'radial' and 'angular' degrees of freedom, respectively. In 
the absence of time-reversal symmetry ('unitary case,' 
(3 = 2) the are arbitrary; whereas the presence 

of time-reversal symmetry ('orthogonal case,' (3 = 1) 
imposes the additional constraint S T = S and hence 
w (3) _ U ( 1 ) T ; v ( 4 ) — ■i/ 2 ) T £j Depending upon the sym- 
metry, the manifold of S-matrices has the dimension 



da := 



2N 2 
47V 2 



N : j3 ■ 



The manifold is endowed with a metric 
g(dS 1 ,dS 2 ) := tr(dSjdS 2 ) 



(9) 



(10) 



which is invariant under the group action dSi — > UdSiV 
(U, V unitary). With respect to this metric, radial vari- 
ations (r a — > r a + dr a ) and angular variations (v^ — > 
v^ 1 ' + dv^) of the S-matrix are orthogonal; the metric 
tensor (g^u) is thus built of two blocks, one 'radial' and 
one 'angular.' The radial part of the metric tensor reads 



9ab ■= g(dS/dT a ,dS/dn) 
" 2r (l-r ) do4 ' 



with inverse 



2r a (l-T a )6 ab 



(11) 



(12) 
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The metric defines an invariant volume element (Haar 
measure) 



dfi(S) = v / [det^|n dx ' 1 

OC 

=^({r})n dr a-n d [ an s ies ] 



where Jg is the Jacobia: 



ni 



J a = 



n„< b k-T 6 |-n c i/v^ -.0 = 1 



Yla 



= 2 



(13) 



(14) 



(The exact form of the angular factor need not concern 
us.) This volume element in turn permits the definition 
of probability densities P(S), with normalization 



dfj,(S)P(S) = 1 



(15) 



In the above parametrization {r a ,v^}, the conduc- 
tance can be easily expressed as a function of S. Accord- 
ing to Landauer's formulaEH it reads 



G= (2e 2 /h)J2r a 



(16) 



(with the factor 2 accounting for spin). The sum T := 
Y^ a r a is called the dimensionless conductance. Assum- 
ing a uniform probability density on the manifold of S- 



matricesJHj P(S) = const, one obtains the well-known 



results 



IT\ P N \ N A N 



and 



var^(T) 



n(n+i) 2 

(2N+l) 2 (2N+3) 

N 2 
4(4JV 2 -1) 



(3=1 
(3 = 2 



(17) 



(18) 



For N — > oo, these imply the universal magnitude of weak 
localization, (T)i — (T) 2 — > (—1/4), and of conductance 
fluctuations, var j3 (T) -> (1/8(3). 



III. BROWNIAN MOTION MODEL 

In the Brownian motion model the probability density 
is no longer uniform and acquires an explicit dependence 
on the parameter X. As X is being varied, X — > X + 
AX, the system executes a random walk on the manifold 
of S'-matrices and, consequently, the probability density 
spreads according to a diffusion equation. In order to find 
the form of this diffusion equation, let us consider - as a 
simple example - a random walk on a flat d-dimensional 
manifold with Cartesian coordinates {x* 1 ,^ = l...d}. 
At each step the velocity v of the moving 'particle' is 



random, with an isotropic distribution so that 
and (v^v") = (l/d)(v 2 )S^. Hence 

P(X + AX) = P(X) + ^(AX) 2 (v 2 ) A P(X) + . . 



(19) 



Upon identifying a fictitious 'time' t := (AA) 2 and the 
diffusion constant D := (l/2d)(v 2 ) , we arrive at the dif- 
fusion equation dP/dt = DAP. This is easily general- 
ized to the (curved) manifold of S-matrices. There, too, 
increments AA are related to a fictitious 'time' t by 



t = (AXf 



(20) 



(Strictly speaking, this relationship holds only for in- 
finitesimal parameter values; the exact relationship be- 
tween t and finite values of the physical parameter re- 
mains unknown.) The diffusion equation then reads 



d_ 
d~t 



P = DpApP 



(21) 



Here Ap is the Laplace-Beltrami operator (the general- 
ization of the Laplace operator on curved Riemannian 
manifolds), and Dp is the diffusion constant 



1 

2d, 



dS dS 
dX' dX 



(22) 



Both are labelled by the symmetry index (3. The aver- 
age (g(. . .)) is a measure for the typical step size of the 
random walk; it in turn defines a characteristic 'time' 
scale 
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dS dS 
dX'dX 



(23) 



In order to make use of this model, one needs to know (i) 
the explicit form of the Laplace-Beltrami operator and 
(ii) how diffusion affects the conductance. 

The Laplace-Beltrami operator on an arbitrary Rie- 
mannian manifold with coordinates {x^} is given bycil 



E 







J V|det, 



| — 



(24) 



Its explicit form on the manifold of S'-matrices, with 
parametrization {r a , v^ 1 '}, follows directly from the prop- 
erties of the metric. As the metric tensor is composed of 
a radial and an angular block, the Laplace-Beltrami op- 
erator can be split into two parts, 



A fJ = A P , T + A 



(25) 



the first containing derivatives with respect to r only, 
and the latter containing derivatives with respect to an- 
gles only. The two parts are called radial and angular, 
respectively. When applied to the conductance, only the 



3 



radial part Ap^ contributes. This radial part is obtained 
by inserting (12) and ( [l3l) into (p4|): all r-independent 
factors cancel, leaving 



^ J/3 OT a 



d 

r Q (l - T a )Jp — 



(26) 



The Laplacc-Beltrami operator has two important prop- 
erties. First, it is Hcrmitian in the sense 



dfi(S) fAph = / d/i(5) hA f 



(27) 



provided / and h are sufficiently smooth. Second, the de- 
viation of the dimensionless conductance from its mean, 
SpT := J2b Tb ~ C0/8> ^ s an eigenf unction: 



2 Et 



<9 



J/3 9r h 



Tb(l - Tfe) J/3 



(3N< 



2^(1 -2r„) 

a 

-{d /N)5 T 



2^r a (l-T a ) 



9 In J/3 



(28) 



Since the Laplacc-Beltrami operator is Hcrmitian, we 
are free to describe the evolution of expectation values 
(/) = / dfiPf either in the 'Schrodinger picture,' P — > 
P(t), or in the 'Heisenberg picture,' / — > /(i). Choosing 
the latter, we see that diffusion causes the conductance 
to approach its equilibrium value exponentially: 



S G(t) := exp[tD A ] SpG 

= exp[-t/(2Nto)]6f}G 



(29) 



This immediately yields the conductance-conductance 
correlation function 



C(t) = (6 G5 G(t)) 

= C(0) • exp[-t/(2iVto)] 



(30) 



The shape of the weak localization peak follows from a 
similar consideration. At t = B = the probability den- 
sity is non-zero only on the submanifold of time-reversal 
invariant S-matrices: P(0) oc S(S T — S). As B and hence 
t increase, this distribution diffuses over the entire man- 
ifold of S'-matrices. We thus find 



W(t) = / d/i(S) P(t) S 2 G 



d(x(S) P(0) 5 2 G(t) 
W{0) ■ exp[-t/(2Nto)] 



(31) 



Equations (|30|) and (^Tj) together imply the key result of 
this paper: 



C(t)/C(0) = W{t)/W(0) = exp[-t/(2Nt )} 



(32) 



IV. DISCUSSION 

We are led to the following three conclusions, (i) The 
conductance-conductance correlation function in B has 
the same shape and width as the weak localization peak: 



C(B)/C(0) = W(B)/W(0) VP 



(33) 



This is a stronger version of the approximate equality 
(||). Which shape the functions have, depends on the 
relationship between t and B. For finite values of P, this 
relationship is not known, (ii) For small P, however, we 
know from equation (E(l) that t ~ P 2 ; hence 



C(P)/C(0) = W(B)/W(0) ~ 1 - C(P 2 ) 



(34) 



This is not a very strong result, but it at least excludes 
linear line shapes, of the kind which are observed in non- 
chaotic cavitiesfl (iii) Provided we are in the regime 
wheije-J ~ P 2 , the width of both functions increases as 

The first conclusion is consistent with semiclassics 
only in the limit P — ► or N — > oo (which implies 
a c \ — > oo); the second conclusion is trivially correct; 
and the third conclusion agrees with Hamiltonian ran- 
dom matrix theorycj except for very small N. It there- 
fore seems that our Brownian motion model is exact in 
the limit of small parameter values and a large number 
of channels. Outside this limit, however, the model can 
describe the phenomena at hand only qualitatively, not 
quantitatively. 

The partial failure of the Brownian motion model re- 
veals that the 'walk' on the manifold of S'-matrices, in- 
duced by varying an external parameter such as P, is 
not really 'random': it cannot be characterized solely by 
its typical step size. In other words, the second moment 
(g(dS/dX, dS/dX)) is either not appropriate or not suf- 
ficient, or both, to characterize parametric motion; there 
must be other constraints O To find the nature and phys- 
ical origin of these constraints, should be the goal of fu- 
ture investigations. 
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